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FINITE ELEMENT METHODS FOR THE STOCHASTIC 
ALLEN-CAHN EQUATION WITH GRADIENT-TYPE 
MULTIPLICATIVE NOISES 

XIAOBING FENG*, YUKUN lO, AND YI ZHANG* 

Abstract. This paper studies finite element approximations of the stochastic Allen-Cahn equa¬ 
tion with gradient-type multiplicative noises that are white in time and correlated in space. The 
sharp interface limit as the parameter e —)■ 0 of the stochastic equation formally approximates a 
stochastic mean curvature fiow which is described by a stochastically perturbed geometric law of the 
deterministic mean curvature fiow. Both the stochastic Allen-Cahn equation and the stochastic mean 
curvature fiow arise from materials science, fiuid mechanics and cell biology applications. Two fully 
discrete finite element methods which are based on different time-stepping strategies for the nonlin¬ 
ear term are proposed. Strong convergence with sharp rates for both fully discrete finite element 
methods is proved. This is done with a crucial help of the Holder continuity in time with respect to 
the spatial L^-norm and if ^-seminorm for the strong solution of the stochastic Allen-Cahn equation, 
which are key technical lemmas proved in paper. It also relies on the fact that high moments of 
the strong solution are bounded in various spatial and temporal norms. Numerical experiments are 
provided to gauge the performance of the proposed fully discrete finite element methods and to study 
the interplay of the geometric evolution and gradient-type noises. 

Key words. Stochastic Allen-Cahn equation, stochastic mean curvature fiow, gradient-type 
multiplicative noises, phase transition, fully discrete finite element methods, error analysis. 
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1. Introduction. The Allen-Cahn equation refers to the following singularly 
perturbed heat equation: 

(1.1) Ut-Au+\f(u)=0 mVT:=Vx{0,T), 

where P C (d = 2, 3) is a bound domain, T > 0 is a fixed constant. Equation 
113 was proposed by Samuel M. Allen and John. W. Cahn in [1] as a mathematical 
model to describe the phase separation process of a binary alloy quenched at a fixed 
temperature. In the equation, u represents the concentration of one species of the 
binary alloy mixture. The positive small parameter e is known as the interaction 
length of the transition region between the domains {x e V : u{x,t) « —1} and 
{x ^ V : u{x^t) ~ 1}, and / is the derivative of a smooth double equal well density 
function F taking its global minimum 0 at rt = ±1. A typical example is 

(1.2) f{u)=u^—u and F {u) = ^ 

which will be used in this paper. We note that t in the equation is a fast time, 
representing ^ in the original Allen-Cahn equation. 

To uniquely determine a solution of the Allen-Cahn equation, a boundary and a 
initial condition must be prescribed. In this paper we shall consider the following set 
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of the boundary and initial conditions: 


(1.3) 7 ^ = 0 in dVr := dV x (0, T), 

on 

(1.4) u = uq inV X {0}, 

where n denotes the unit outward normal to dV and ^ = \/u • n denotes the normal 
derivative of u. 

Besides the important role it plays in materials science, the Allen-Cahn equation 
has also be well known because of its connection to a well-known geometric moving 
interface problem called the mean curvature flow, which is governed by the following 
geometric law: 

(1.5) Vn{x,t) =-H{x,t) xeTf, 

where and respectively stand for the normal velocity and the mean 

curvature of at the space-time point (x,t). It was proved in [20l[6l[T5] that when 
e ^ 0, the zero-level set := {x £ V : u{x,t) = 0} of the solution of problem 
(fnt-drit converges to the mean curvature flow {rt}t>o- Furthermore, because of 
the above property the Allen-Cahn equation has become a fundamental equation and 
a building block in the phase field methodology for general moving interface problems 
Eilll. We refer to [H |T6] and the references therein for an account on numerical 
methods for problem 

For application problems of the mean curvature flow, there may exist uncertainty 
which arises from various sources such as thermal fluctuation, impurities of the ma¬ 
terials and the intrinsic instabilities of the deterministic evolutions. Therefore, it is 
interesting and necessary to consider the stochastic effects, and to study the impact 
of the noises on regularities of the solutions and their longtime behaviors. This leads 
to considering the following stochastically perturbed mean curvature flow [131 Ell ES] : 

(1.6) Vn{x,t) = -H{x,t) ^Sx{x,t) • n x eTt, 


o 

where X : x [0, T] ^ and X denotes the Stratonovich derivative of X, (5 > 0 is 

the noise intensity that controls the strength of the noise. Formally, the phase field 
formulation of ( |1.6[ ) is given by the following stochastic Allen-Cahn equation [21]: 

fill 1 O 

(1.7) — = Au-^f{u)+5Wu-X, 

at 

More specifically, we assume X is a vector field-valued Brownian motion that is white 
in time and smooth in space, i.e.. 


( 1 . 8 ) 


X{x,t) =X{x)W{t), 


where X : is a time independent deterministic smooth vector field with 

compact support in V and W(t) is a standard M-valued Brownian motion (i.e., Wiener 
process) on a given filtered probability space : t > 0},P). Thus the 

Stratonovich stochastic partial differential equation (SPDE) (ini) becomes 


(1.9) 


du = 


Au-\f{u) dt + 6Vu-XodW{t), 
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and the corresponding ltd SPDE is given by 


( 1 . 10 ) 


du = 


1 (5^ 1 

Am - -^f{u) + — V(Vm ■X)-X dt + 6Vu ■ X dW{t) 

C Zi 

1 s‘^ 1 

+ —{B: D‘^u + 6 • Vrt) dt + 5Vu • X dW{t) 


1 


2 


-^/('^) H-(div(5Vi4) + (6 — div5) • Vi^) 

6 .Z 

^5Vu- XdW(t), 


dt 


with hj = VXj ' X 


where B = X^X e with Bij = XiXj, b = (VX)X e 
and div5 — b = (divX)X. Note that the ltd SPDE ( |1.1Q[ ) has two extra terms which 
are absorbed in the Stratonovich SPDE (1.9). The stochastic Allen-Cahn equation 
(1.10) will be complemented by the boundary and initial conditions (1.3) and ( |1.4[ ). 

Remark 1. We can also consider the following more general vector field /^T] / 

m 

%{x,t) = Y,^\x)Wk{t), 

k=l 

where m is a positive integer, Wk{t) are independent R-valued standard Brownian 
motions, and X^ : ^ is a time independent deterministic smooth vector field 

with compact support in V for k = 1,2, • - ,m. In this case, the corresponding B and 
b in the SPDE (1.10) are given by 

m 

B 




and 


b = Y^{VX’^)X'^ 


k=l 


k=l 


We remark that the results of this paper can be easily extended to the case with the 
above general X. 

In this paper we assume that 


^0 e C^{V), 

'^4,/3o ( 


( 1 . 11 ) 

(1.12) A e Co’'’“(D) 

for some G (0,1]. It was proved in [HJ Theorem 4.1] that if the domain V is 

hold, there exists a unique strong solution 


smooth and the assumptions (1.11)-(1.12 
u such that 


(1.13) 


(1.14) 


nt , 

a(x, t) = rto(^) + / Au{x,s) - ^f{u{x,s))-^—(B:D‘^u{x,s) 

Jo ^ ^ ^ ^ 

-\-b ■Xu{x,s)^ dsJ 6Xu{x,s) ' X{x) dW{s), 


?=» 

on 


on OVt 


hold P-almost surely. Moreover, u{',t) is a continuous C^’^(P)-semimartingale for 
any D < p < po- Eurthermore, for any multi-index |cr| < 3 and p > 1, there exists a 
positive constant Cq = C{p, 6, e) such that the p-th moments of the spatial derivatives 
of u satisfy 


(1.15) 


sup E 


te[o,T] ^xev 


sup \d^u{x, t)|^ 


<Co, 
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where and throughout this paper E['] denotes the expectation operator, 


We note that other solution concepts for the SPED problem (1.10), (1.3)-(1.4) 


were also studied in the literature, this includes mild solutions based on the semigroup 
approach [18] and variational solutions El based on the variational approach. In this 
paper we shall assume that there exists a unique strong solution to problem (1.10), 
(1.3)-(1.4) satisfying (1.13)-(1.15). It is clear that the strong solution is also 

an adapted i7^(P)-valued process such that for any t G (0,T] there holds P-almost 
surely 


(1.16) 


4>) = («o, <P)- [{I + yS) V«(s), V0) ds 

- 4 /* ifiHs)), <^) ds - 4 f ((divB - b) ■ Vw(s), (/.) ds 

+ 6 [ {Vu{s)-X,(f>)dW{s) 

Jo 

where (•, •) denotes the standard inner product on L‘^{V). 

In addition to the results about the well-posedness of strong solutions and the 
regularity estimates given by (1.15), the tightness of the solution in the sharp interface 
limit as e ^ 0 was also established in [21] for the stochastic Allen-Cahn problem 
ra , ([n -( |1.4 ). However, the rigorous justification of the convergence of the zero 
level set of the solution u to the stochastic mean curvature flow (1.6) is still missing. 


In [ain], the stochastic Allen-Cahn equation in one dimension with additive space- 
time white noise was studied and the sharp interface limit was established. Finite 
element methods for this one dimensional model was also proposed in [12]. But in 
higher dimensions, it requires spatial correlations for the noise otherwise the space- 
time white noise is too rough to allow the existence of solutions in reasonable function 
spaces. In m, finite element methods for the stochastic mean curvature flow of planar 
curves of graphs were investigated, where the SPDE is quasilinear and arises from the 
level set formulation of the mean curvature flow. However, it is not clear if the results 
of [To] can be extended to higher dimensions. Very recently, finite element methods 
for a stochastic Allen-Cahn equation with Lipschitz continuous multiplicative noises 
were proposed in m, the strong convergence with suboptimal rates was proved. 

The primary goal of this paper is to develop and analyze two fully discrete finite 
element methods for approximating the strong solution of the stochastic Allen-Cahn 


problem (1.10), (1.3)-(1.4) and to establish strong convergence with optimal rates 


for the proposed fully discrete finite element methods. Our strong convergence result 
will be established under the assumption that the strong solution G 

{s > 3) for any t G [0,T] and it satisfies the following estimates: 


(1.17) 


sup E 
te[o,T] 


|u|r 






<Co = C{p,S,e) Vp>l. 


We note that the above assumption is reasonable in view of (1.15). Since we are 


interested in the case where e is small, we shall make an effort to track the dependence 
of all constants on the parameter e in our convergence analysis. We remark that for 
the deterministic Allen-Cahn equation, using a nonstandard technique it was proved 
that the error estimates depend on ^ only in some lower polynomial order mm- 
Such error estimates can be used to prove convergence of the numerical interface to 
the mean curvature flow. However, it is not clear whether the techniques used in 
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[m [3 |9] can be extended to the stochastic case. As a result, our error estimates 
are shown to depend on ^ exponentially. Another goal of this paper is to use the 
stochastic Allen-Cahn equation as a testbed to develop numerical analysis techniques 
and machineries which hopefully are applicable to other stochastic PDEs. To the best 
of our knowledge, there is no numerical analysis result for nonlinear stochastic second 
order PDEs with gradient-type multiplicative noises in the literature, this is partially 
because the gradient-type noises in a SPDE like (1.10) give new challenges for error 


estimates. More specifically, such a multiplicative noise introduces additional diffusion 
and convection terms to the PDE, and the resulting SPDE could be convection- 
dominated if the multiplicative noise is too strong. To avoid the latter situation, we 


shall assume that S is not too large in (1.10). 


The remaining of the paper is organized as follows. We present in Section some 
technical lemmas which play a crucial role in our convergence analysis. In Section 
we propose two fully discrete finite element methods for the stochastic Allen-Cahn 
problem ( 1.10[ ), ( |1.3[ )-( pr4| ). These two methods differ only on how the nonlinear term 
is discretized in time, which leads to different conditions on the time step r for the 
well-posedness and error estimates. We establish the strong convergence with optimal 
rates for both fully discrete finite element methods. Einally, in Section we present 
some numerical experiments to illustrate the performance of the proposed numerical 
methods and to verify the theoretical error estimates obtained in Section as well as 
to examine the dynamics of the numerical interfaces, in particular, pertaining to the 
noise strength. 

2. Preliminary results. Standard space notation will be used in this paper [4], 
in particular, we point out that H^{V) for /c > 0 denotes Sobolev space of order k 
and (•, •) denotes the standard inner product of L‘^{V). Throughout the paper C will 
be used to denote a generic positive constant independent of the parameters e, 6, and 
the space and time mesh sizes h and r, which can take different values at different 
occurrences. 

The goal of this section is to establish several technical lemmas that will be 
crucially used in Section Each of these lemmas shows a certain property of the 
strong solution u of the stochastic Allen-Cahn problem ( 1.1Q[ ), ( |1.3[ )-( pr4| ), such as 
(e, (5)-explicit energy estimates and Holder continuity in time with respect to the 
L‘^(V)- and i^^(P)-norm. 

We begin with the following uniform estimate for the expectation of the p-th 
moment of the Cahn-Hilliard energy functional 


Av) := dx. 


A similar result can be found in [21] without tracking the dependence on parameters 
in the estimate. We include a proof here for completeness. 


Lemma 2.1. Let u he the strong solution to problem (1.10), (1.3)-(1*4). There 
holds for any p>l 


r 

(2.1) sup E[j(u(i))^]+E / pJ{u{s)r~^\\w{s)\\i^, .ds 

te[o,Ti L^o ^ ^ 


<Ci, 


where 


w := —Au + 


and 


Ci=e 


- 

















Proof. By direct calculations we have 
D{J{u)P) = (-Am + ^F' 

D\j{uf) = p{p - l)J{uf-^ (-Am + 0 (-Am + ^F' 

+ pJ(w)P-i(-A+^F"(u)). 

Now applying Ito’s formula to the functional := J{u{')Y and integration by 

parts, we obtain 

(2.2) J{u{t)Y = J{uoY - ( pJ{u{s)Y~^\\w{s)\W 2 (T::)ds 

Jo 

+ pJ{u{s)Y~^ (m^(s), y {B ; D‘^u{s) + h ■ Vu(s))) ds + Mt 

+ y ^ P{P- l)J{u{s)Y~'^{w{s),Vu{s) ■ Xf ds 

f pJ{u{s)Y-^ [ [|V(V«(s)-A)|2 + 4f"(«(s))(V«(s)-A)' 

Jo Jv^ ^ 


S2 rt 

where Mt is the martingale given by 


dx ds, 


(2.3) Mt=S [ pJ{u{s)Y~'^{wis),Xu{s)-X)dW{s). 

Jo 


By integration by parts and a direct calculation, we have 




(2.4) y [ [|V(Vti(s) • X)|2 + 4F"(ys))(Vw(s) • A)' 
Jv ^ 

(52 


dx 


2 jt> 

+ (■*^( 5)1 y {B ■ D'^u{s) + b ■ Vit(s)) j 
= y y J (g{x) : (Vit(s) 0 Vit(s)) + ^g{x)F{u{s))j dxds, 

where G{-) : V —y anci g{-) : V —> M are (iefineci by 

(2.5) Gij = [dk{XkdiXi)]5ij + di[XkdkXj] - 2dk[XkdiXj] + {dkXiYdkXY, 

(2.6) g = dkiXkdiXi], 


Taking expectation on both sides of (2.2), by (2.3) and the fact that E[M(] = 0, 
we get the following equation 

r rt 

(2.7) E [J(M(t))P] = J(uo)P - E 




y pJ{u{s)Y ^ iy(s)|li 2 (p) ds 
y pJ{u{s)Y~^ y (g{x) ; (Vn(s) 0 Vn(s)) 


1 

+ ^ 


dx ds 


S 2 r .t 

+ -E 


/ p(p-I)J(u(s))'’-^(w(s),Vu(s)-Xf ds . 
Jo 


6 














Now it remains to estimate the third and fourth terms on the right-hand side of 


(2.7). On noting that 

the third term can be bounded as follows 


52 

(2.8) -E 


pJ(u(s))P-i J (g{x) ; (Vu(s) 0 Vu(s)) -h ^g{x)F{u{s))'j dxds 

<Cp5^X\\l,^^^ l\[J{u{s)r] ds. 

For the fourth term, by integration by parts and the fact that 


we have 
(2.9) 


^2 r .t 


f p{p — l)J{u{s)y {w{s),\/u{s) ' X)‘^ ds 
Jo 

< Cp{p-l)6^X\\l,^^^ J\[j{u{s)r] ds. 

Finally, combine (2.7)-(2.9) and apply the Gronwall’s inequality, we obtain the 
desired estimate (2.1). □ 

Remark 2. In the case ofp > 1 in Lemma\2^ the fourth term on the right-hand 
side of (2.7) can also be hounded by 

£2 r rt 

(2.10) —E / p{p-l)J{u{s)Y-^{w{s),Vu{s)-Xf ds 

2 L^o 

<{p-l)5^\\X\\l^^^E J\j{u{s)r-Y\w{s)\\l,^j,^ ds . 


Therefore (2.1) is replaced by 

( 2 . 11 ) 


sup E [J(u(f))P]-h OpE 
te[o,T] 


[ pJ{u{s)Y ^ ||w(s 
.Jo 


Il2(d) 


ds 


<Gi, 


where 

Cp = l-{p- l)5^\\X\\l^^^ and Oi = [J(uo)]'’. 

However, ( |2.11[ ) implicitly imposes a restriction on 6 with respect to p, namely, 5 < 

1 /{Jp ^\^\\c{f)))' 

Remark 3. Note that in the case ofp = 1 and 6 = 0, we recover the the following 
deterministic energy law: 

sup J{u{t))F / \\w{s)\\^ 2 rj:)^ ds < J{uo). 
te[o,T] Jo 
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Next we derive Holder continuity in time for the solution u with respect to the 
spatial I/^-norm and the H^-semmoim. Both results play a key role in our error anal¬ 
ysis in Section In a sense they substitute the time derivative of the solution, which 
does not exist in the stochastic case, appeared in the deterministic error analysis. We 
also note that, as expected, the constants in these estimates depend on in some 
polynomial orders. 


Lemma 2.2. Let u he the strong solution to problem ( 1.1Q[ ), ( |1.3[ )-( pr4| ). For any 
5 , t G [0, T] with s < t and p>2, we have 


( 2 . 12 ) 




r ik(C) - «wiiii(%iiv(t^(c) - dc 

J s 

f ||«(C) - «(s)|liiy)||v(«(C) - «(s)) ■ dc 

J s 

||u(t) <C 2 {t-s) 


TT. 

— E 
4 


+ E 


for some C 2 = C{p,5, \\X\\c 2 ^^^,T)€-p/^{[J{uo)]p/^ + C) > 0. 


Proof. Applying Ito’s formula to the functional ^ (?/(•)) := \\u{')—u with 

fixed s G [0, T) and using integration by parts, we obtain 


(2.13) \\uit)-uis)\\l,^^^ 

=pf^ iit^(c) - u{s)rr.i^^dc 

+ ^ ^ l|w(C) - M(s)lli 2 (% {u{0 - u{s), B : D‘^u{C,) + b ■ V«(C)) df 

+ pf llt^(C) - t^(s)lliA) («(C) - u{s),SXuiO ■ x) dWiC) 

+ ^ jfp~ 2)ll«(C) - «(s)) («(C) - u{s), Vu{C,)-xf dC 

+ ? / 11^(0 - «(s)|liiy)||Vw(C) • dC 

= I + II + III + IV + V. 


We only consider the case p > 2 since the proof for the case p = 2 is simpler. By 
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Young’s inequality, we have 


(2.14) E 


U(C) -u{s 


-E 

-E 


< -^E 
- 2 

+ -E 


< -^E 
- 2 


pj'^ h(C) - ti(s)||i7y)||Vti(C) - Vw(s)||i.(^) dC 
pj'^ ||«(C) - «(s)|liiy) (Vti(C) - V«{s), Vti(s)) dC 
f llt^(C) - uis)\\ll^^^\\Vu{0 - Vuis)\\l,^r>) dC 

.J S 

.J s 

||«(C) -«(s)|li;y)||v«(C) - v«(s)||i.(^)dC 
" /*E[|KC)-«(s)||i.(J dC + E 

J S ^ 


p — 2 


l|VM(s)||i2(i,) (t-s) 


Similarly, we have 


(2.15) E 


pj^ ll«(C) - («(C) - «(s), -^fHO)) dC 


WT^v) 


< 


[ HO-4 

J s 

X (^11404(0 - u{s))\\h(v) + ^l4(w(0)llLh®)) dC 
2(p-2)^ E II4C) -40lli2(o) 


dc 


+ Cpe-P/^( sup E[||u(C)||f4(c)l + sup e[ 4h40) 


s<C<t 


p/2 

Li(7^)J 


Using (2.14)-(2.15), the first term on the right-hand side of (2.13) can be bounded 
as follows 


(2.16) E[/] < -|E 


f 4(0 - u(s)4l(%||Vu(0 - Vu(s) 114(1,) dC 

J s 

54-2) [' 

2 Js 


E 


4(0 - w(0II4(d) 


dC + E 


l|Vu(s)||4(i,) (t-s) 


1 


+ C,e-P/^( sup E [4(0114(^)1 + sup E ^F{u{C)) 
^s<C<t ^ s<C<t ^ 


p/2 

LiCD)J 




Again, by using Young’s inequality and integration by parts the second term on 
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the right-hand side of ( |2.13 ) 
5^p 

(2.17) E[//] < — 


can be bounded by 


^V(C)-t^WlliI(P)l|Vn(C) 


•X- 




+ 


o 

YiWXfciV) + [l|V«(C)lli 2 (i,)] (t-s) 


Next, it follows from Young’s inequality that 

^yP(p - 1) f 11^(0 - «WlliI(%l|Vt^(C) • dC 

- l)(p - 2) E [mo - dC 

4- ^r. _ 1U2|| Flip _ OUT, IFT rilTyi/r^MIP^ 


(2.18) E[IV + V]<E 
< 


+ {p- 1)(5^||X||^(^^ ^sup^E [||Vn(C)||i.(^)] {t - s). 

Combining (2.13), (2.16)-(2.18) and using the facts that < Se^J{u) -\-C 

and E[III] = 0, we have 


(2.19) 


E 


+-E 
2 

5^p 


< 


\u{t) - u{s)\\l2^x>~) 

f\\ui0-uis)M^^^\\V{u{0-u 

.j S 

E [ f ||«(C) - «(s)||^I(%||V(«(C) - U{s)) ■ X\\l 

.J s 

2 '-{pd‘^ + 0 E lk(C)-«(s)lli2(i,) dC 
C{pO\\\XM^vM-p'H sup E [j(«(C)^’/2)l +c)(t 


P-2, 


iO-u{s))-X\\l.^^^ dC 

dC 

-s). 


Finally, the desired estimate (2.12) follows from (2.19), the Gronwall’s inequality 
and Lemma |2]^ The proof is complete. □ 


Lemma 2.3. Let u he the strong solution to problem (1.10), (1.3)-(1*4). For any 
s,t E [0,T] with s <t, we have 


(2.20) E[\\V{u{t) - «(s))||i 2 (p)] + ^E ||A(«(C) - u{s))M^^) dQ 


< Cz{t- s), 


where 


C3 = C S‘ 




+ ;2 + ^ 


Cl + sup E 
s<C<^ 


llAn(c)iii2(p)]). 


Proof. Applying Itb’s formula to the functional ^(i4(-)) := ||Vi4(') — Vi4(s) ||^ 2 ('x)) 
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with fixed s G [0, T) and using integration by parts, we get 

(2.21) \\^u{t) - Vm(s)|||2(x)) = (^■“(0 - Am(s), Au{C)) dC, 

+ 2 1* (Aw{C) - Ati(s), ^/(w(C))) dC 
- 6^ (a^/(C) - A^/(s), B : DMO + b ■ Vm(C)) dC 

-26 f (Au(C)-Au(s),Vu(()-X)dW(C) 

J S 

\V{Vu{C)-X)\^dxdC. 

The first term on the right-hand side of ( 2 . 21 ) can be estimated by Schwarz 
inequality as follows 


( 2 . 22 ) 


2E 


L’y s 

r rt 

= -2E 


f (Au(C) - Au{s), Au{C)) dC 

J S 

[ I|Am(C) - A'u(s)||| 2 (i,) ciC+ / {Au{C) - Au{s),Au{s))dC 

J s J s 


< -E 


LJ s 
t 


[ I|Am(C) - AM(s)|||2(i,)ciC 

y s 


l|AM(s)||i2(i,)(i-s). 


For the second term on the right-hand side of ( 2 . 21 ), by Sobolev embedding ^ 

L^{V) for d < 3, we have 


(2.23) 


2E 


(Am(C) - A'u(s), ^/(m(C))) dC 

^ (||Am(C) - Am(s)|||2(i,) + ||/(«(C))lli2(i,)) dC 


2 ^ 

< 


<^( sup E 


e ^s<C<t 


l|Au(C)||i2(^) 


+ 


sup E[J{uiC)f]+c){t-s). 


s<C<t 


Next we bound the third and fifth terms on the right-hand side of ( 2 . 21 ) as follows 

rt 


(2.24) 


E 


52 ^ (Au(C) - Au(s), B : D'^u{0 + b • Vu(C)) dC 

J ||Au(C) - AM(s)|||2(i,)(iC +^<^'‘11^11^(5) 

X ( sup E [||Au(C)||i 2 (i,)] + sup E [||Vu(C)||i 2 (i,)])(i-s), 


< -E 
“ 2 


s<C<t 


and 

(2.25) 


E 


\X{Xu{0-X)\^dxdC^ <C5^\\X\\%^^^ 


X ( sup E 


ll^'^(C)lli2(7:») 


■ sup E 
s<C<t 


l|Vu(C)lli2(c)])(t-s). 
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The desired estimate (2.20) then follows from (2.21)-(2.25), Lemma 2.1 and the fact 
that the expectation of the fourth term on the right-hand side of (2.22) is 0. □ 

Finally, we include a Holder continuity estimate for the nonlinear increment 
f{u{t)) — f{u{s) in the L^-norm which will be useful to control the nonlinear term in 
the error analysis. 

Lemma 2.4. Let u he the strong solution to problem (1.10), (1.3)-(1*4). For any 
G [0,T] with s <t, we have 

(2.26) E[||/(«(i)) - /(m(s)) 111,2(1,)] < Ciit - s), 
where C 4 = C{S,e, \\X\\ci(t>)^Co,Ci). 

Proof. Applying Ito’s formula to the functional 4 >(i 4(')) := Ili 2 (x)) 

with fixed s G [0,T), we obtain 

(2.27) ||/(«(i))-/(«(s)) 111 , 2 ( 1 ,) =2 /" f {f{u{C))-f{u{s)))f{u{C)) 

J s 

X Aw(C) - + ydiv(BVw(C)) + y (& - divB) • Vw(C)] dxdQ 

+ 26 f [ {fiuiC))-f{uis)))nu{C))Xu{0-XdxdW{C) 

Js Jv 

(/(■w(C)) - /(«(«)))/"(■w(C))|Vm(C) • X\'^dxdC 

+ r / [f{u{0)f\Xu{0-XfdxdC 

Js Jv 

Taking the expectation on both sides, it follows from integration by parts. Young’s 
inequality and Sobolev embedding ^ L^(V) for d < 3 that 


(2.28) 


E 


Wfiuit)) - /(t^(s))iii2(i,)J < c(i + (5^iix|iy(^))(t-5) 

X rC+(l + e“2) sup E [||■u(C)||ioo(I,)] + sup E [J(m(C))^]). 
^ s<C<t ^ / 


s<C<t 


Finally, the desired estimate (2.26) follows from (1.17), (2.28) and Lemma 


2.1 


3. Finite element methods. In this section we propose two fully discrete finite 
element methods for problem ( 1.1Q[ ), ( |1.3[ )-( p~T| ) and establish their strong convergence 
with optimal rates. 

Let tn = nr (n = 0,1,..., N) be a uniform partition of [0, T] with r = T/N and 
7^ be a quasi-uniform triangulation of V. We consider the finite element space 

V;::={vheH\V):Vh\KeVr{K) MKeTh}, 

where Vr{K) denotes the space of all polynomials of degree not exceeding a given 
positive integer r on K G Th- Our fully discrete finite element methods for problem 
( fn^ , (fot-di^ are defined by seeking an -adapted H^-valued process {u'^} 
(n = 0,1,..., N) such that P-almost surely 

(3.1) (ul+\vn) + T[{l+^-^B)Vul+\Vvh) +T^ir+\vn) 

= {ul, Vh) - Tj ((divB - 6) • Vul Vh) + 6 (Vu^ ■ X, v^) AW^+i 
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for all Vh € Vf[, where AWn+i ■= W{tn+i) ~ W{tn) ~ A/’(0 ,t) and 


(3.2) 


/"+! := («n+l)3_«n+l 


P+l .= (^n+l)3_. 


We remark that in the deterministic case (i.e., (5 = 0), ( |3.2[ ) corresponds to a fully 
implicit scheme or a convex-splitting scheme. We choose = PhU^ to complement 

(3.1) where Ph : L‘^{P) —is the L^-projection operator defined by 

{PhW,Vh) = {w,Vh) Vh e Vh- 

The following error estimate results are well-known [ill]: 

(3.3) \\w - Phwh^^j,^ + h||V(n; - Phw)\\mv) < 

(3.4) \\w - PhW\\L^(j)-) <Ch‘^~^\\w\\H 2 (T') 
for all w G H^{V). 

We consider a convex decomposition F{v) = F^{v) — F~{v) where 


(3.5) 

and define functionals 


F+{v) = -{v^ + l) 


and 


F-{v) = -v^ 


1 

(3.6) I{v) := -{v,v) + ^i\7v,Vv) + ^i\7v X,\7v X) + ^iF{v),l) 

- {<, ^ ((divi? - b) ■ V<, v) - y (V< • X)AWn+uv ), 

1 q- T(5^ 

(3.7) H{v) -.= -{v, v) + -(Vn, Vn) + — (Vn • X, Vn • X) + ^(F+(n), 1) 

- (l + ^) «, v) + — ((divB - h) ■ V<, v) - (5((V< • X)AW^„+a, ^;). 

It is clear that H{v) is strictly convex for all h, r > 0 and I{v) is strictly convex when 
r < e^. Then a straightforward calculation shows that the discrete problem (3.1) is 
equivalent to the following finite-dimensional convex minimization problems: 


(3.8) 

= a.rgmml{vh), 

when 

/"+! = - ul+^ and T < 




(3.9) 

= argminJJ(v^), 

when 



VHdV- 




Therefore, the existence and uniqueness of problem (|3.1|) is guaranteed for all h, r > 0 
when 

Remark 4. We can also consider a modified scheme 


u'jI and for all r < when 


(3.10) (ul+\vh) + T[{l+^-^B)Xul+\Vvh) +T^{r+\vh) 

+ T^(^{diYB -b) -Xu^+^Vh) 

= {ul,Vh) + S{Vul ■ X,Vh) Aiy„+1 yvh € Vh, 

where we replace the term r^((divR — b) • Vu^^Vh) in ( |3.1[ ) by r^((divR — b) • 
^Vh) in (3.10). Clearly, ( |3.1[ ) has a simpler form and the stiffness matrix is 
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symmetric. On the other hand, (3.10) has one more convection contribution making 
the corresponding stiffness matrix non-symmetric. 

Remark 5. Due to the time discretization, it is unclear whether the discrete 
analog of energy bound in Lemma \2.1\ is valid. However, the error estimate below 
does not require any discrete solution estimate. 

Let e^ := u{tn) — (n = 0 , 1 , 2 ,A^), we now derive error estimates for e^. 

Theorem 3.1. Let u and denote respectively the solutions of problem 

(1.10); (1.3)-( [L4| ) and scheme (3.1); and u satisfies (1.17). In the case of f^^^ = 
~ under the following mesh constraint 


(3.11) 
there holds 


T<C{e-^ + 5^) \ 


(3.12) 


0<n<Ar 


< T 





+ E 


N 

^r||Ve' 

,n=l 

Ci 


,n II2 

LH-D) 


^C(S* + ^)T 


(^3(1 + (5^ +5"^) + 

+ Co(l + 5^ + ^ ^ ^2 min (r+l,s) 

< ^((5, e, T, ||X||, Co, C2, C3, C4) [r + h? . 

In the case of f^~^^ = — uf^, under the following weaker mesh constraint 

(3.13) r < 

we have 


(3.14) sup E 

0<n<N 





+ E 


N 


^r||Ve"||2 


< j,gC(l+6^ + ^)T 


(^?3(1 


,n=l 

<52 + ^4) 


L^V) 

C 2 + C 4 \ 


c -r 


+ Co (1 + <52 + + (Co + 

< 0(5, €, T, ||X||, Co, C2, C3, C4) [r + (r,*-!)] _ 

Proof. We write e'^ = u{tn) — = r]^ where 

:= u{tn) - Phu{tn) and C •= Phu{tn) - ^ 2,..., N. 

It follows from ( |1.16[ ) that for all tn {n > 0) there holds P-almost surely 

rtn+i 

(3.15) {u{tn+i),Vh)-{u{tn),Vh)+ {Vu{s),'^Vh) ds 

,t„ 


+ ^ {fiHs)),Vh)ds+^ {BVu{s),Wvh) ds 

^2 rtn+i 

+ y / ((divB - b) ■ Vu{s),Vh) ds 

J tn 

5 f''^' {Vu{s) ■ X, Vh) dW{s) Wvhe Vh. 

Jtn 
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Subtracting (3.1) from (3.15) and using the decomposition of we obtain the 

following error equation: 


(3.16) (^"+1 -{71^+^ - v^,Vh) - -Vul+\Vvh) ds 

~\f {f(.Hs)) - r^^,Vh)ds 

^ Jtr, 

r2 ptn-\-i 

- - / (B(Vm(s) - V<+i), Vvh) ds 

J tn 

- y ((divB - b) ■ (Vw(s) - VO, Vh) ds 

+ 6 {{Vu{s) - VO • X, Vh) dW{s) yvhG v;;. 


Setting Vh = in (3.16), we get P-almost surely 


(3.17) (r+^ - r,r+^) = ^ (vw(s) - vot ds 

1 rtn+i 

{f{u{s))-r+\e+^)ds 

^ -Itn 

- - / {B{Vuis) - VO'). VO') ds 

Jtn 

_ ^ ((divB - b) ■ {Wu{s) - VO, O') ds 
+ S ((V«(s) - VO ■ X, O') dW{s), 

;= Ti + T2 + Tg + T4 + Ts + Te- 


By the elementary identity (a — b)a = |(a^ — 6^) + |(a — b)^, we get 


(3.18) 


E [(O' 


r,0')] 






llO'lli.(i,)-IIO 

we^^-cwhm 


LHV) 


Next, we bound the right-hand side of (3.17). First, since Ph is the L^-projection 
operator, we have E [Ti] = 0. For the second term on the right-hand side of (3.17), 
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we have by Lemma [2^3] and Young’s inequality that 

■ rt-n 


(3.19) E [T 2 ] = -E 
-E 

< CE 


^ (Vw(s) - Wu{tn+i), vr+^) ds 

rtn+1 

Jt„. 




rtn+1 

/ I|Vm(s) - VM(i„+i)|||2(i,) ds 


—E 
16 


Iivr+'iii2(i,) 


-t E 

J|Vr?”+ii 2 ( 2 ,)’ 

r + ^E 

4 

livr+ii 2 ( 2 ,)^ 

-E 

J|vr+'iii 2 ( 2 ,)’ 

T 



T- —E 
16 


iivr+'iii2(p) 


<C'3r2+E[||Vr,"+i||i2(,,)_ 

In order to estimate the third term on the right-hand side of ( |3.17| |, we write 

(3.20) -(/(n(s)) - r+\r+') = -if{u{s)) - 

- (/(m(Wi) - fiPhUitn+l)),^^^) 

-(/(p„n(t„+i))-r+i,r+')- 

By Lemma [ 2 .4[ we obtain 


(3.21) - E [{f{u{s)) - /(n(t„+i)),r+')] 

1 


< -E 

2 Ct 


-h LiE 
2 




2 Ci 


\\f{u{s)) - f{u{tn+l))\\l 2 ^n) 

Next, by ( |1.17| ) and ( |3.3| )-( [3^ , we have 
(3.22) -E [(/(^(Wi) - /(P/i^(tn+i)), r+')] 

x|K+'lli2(J+CtE 


< -E 

- 4Ct 


< 


< 


< 


||'u(t72+l) “h '^(^n+l)-f/i'^(^n+l) P -P/i'^(^n+l) 

X (e [||r?”+ii2(p)])Nc'tE [||r+'1112(2,)’ 

(e ||P/irt(t„-)-i)||^oo(2,) “I” II'*^(^«-i-i)1Il°°(i>) P \P\ ) 


^ i|lioo(2,) 




40- 


X (e [ii??"+'iii2(2,)]) YctE [iir+ii2(i,) 

^ (e [Wv^+^Un)])" +C^E [||r+' 1112 ( 2 ,) 
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When the last term on the right-hand side of (3.20) can 

be bounded by the monotonicity property 

(3.23) -E [(/(P;,^(t^+i)) - r+\r+')] < E [iir+'iii2(^)' 

When ~ term on the right-hand side of (3.20) resu lts 

in an extra term If can be controlled by using Lemma 

follows 


2.2 


as 


(3.24) E [«+i - ui, r+')] = -E [((r/"+' -»?") + (r+' - D, r+')] 

+ E [{u{tn+i) - u{tn),C~^^)] 

1 " [||r/"+i +CtE [||r+ii2(p; 


<-E 

2 Ct 




E 


||^n+l||2 


LHV) 


H-E 

4Ct 


wcwhiv) +c'tE iir+'iii^(i,) 


Combining (3.20)-(3.24), we choose = 1 when to obtain 

(3.25) E in] < + 1e [||r+' Whiv)] ^ ^ WUt^)] ) " 

and choose C\ = when /"“t^ = ~ to obtain 


(3.26) E [Ta] < 


C 2 + C 4 2 


+E 




l^cd) 


t + Ie 


iiriii.(,,) 


I A 


(e [ll^"+'rn.(i,)])^+ ^E [||r,"+i||i.(^) + ||,?"|| 


^n ||2 

Il 2 (D) 


T. 


Similar to the estimate of T 2 , the fourth and fifth terms on the right-hand side of 


(3.17) can also be bounded by Lemma 2.3 and Young’s inequality as follows 
r 2 r rtn-\-i 

(3.27) E[T 4 ] = -yE 


(^2 


/•hn + l 

/ ((Vit(s) - Vit(t„+i) ■ X, • X) ds 

Jtr, 

/ ((Vr?"+^ + • X) ds 

Jtr,. 


< j^iixiiy^yay + I livr+tli^(p) 


16 


+ (5^||X||y^)E 


r--E 


iivr+'-^11^(1,) 


r, 


and 

(3.28) 


<52 

E^s] = - yE 


r 

—E 


r^n+l , . 

/ ^(divP — b) • (Vu(s) — ds 

.Jtn 

^ ((divB - h) ■ (Vr/" + vr), r+^) ds 


< iix|iy+ (1 + iix|iy(25))<5 "e iir+'iii^d,) 


+ llx|iy(g)E ||V7?"iii.(^) 


r H-E 

16 


iivriii.(,,) 
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By the martingale property, Ito’s isometry and Lemma [231 we have 


(3.29) E[T 6 ]<^E ir+'-rili 2 (i,) 


+ —E 
2 




< -E 
“ 2 


ll^n+l _ ^n||2 


+ y(l + C 0 E 


+ —E 
2 


L 2 (D) 

rtn-\- 


(52 


< ^E 

- 2 


||(Vr," + Vn-^lli=(^) 

(52 


T H-E 

2 C' 


ii(V7?” + vr)-^iii.(i,) 


iir+' - cwUv) + y(i + c')\\x\\i^^^CsT^ 


+ ^)\\X\\%^.E 


2 

\C(vy 




+ yE [||vr • X|| 2 ,(^)] r + 6^(^ + ^] ||X||^(g)E [\\VC\\hm_ 


Now taking C' = and C" = in ( |3.29[ ), we obtain an estimate 


for the last term on the right-hand side of (3.17) 


cm 


(3.30) E[Te] < -E ||^"+i -^"i|2 


L^V) 


+ yd + 16i"||.Y||J,,g,)||.Y|£,g,(73r" 




+ 


16^2||X||2 


cm 


)ll^llc(S)E 




(52^ 


iivr-^ 111 .( 1 ,) 


"+I 6 ® 


mmv) 


Taking expectation on (3.17) and combining estimates (3.18)-(3.19), (3.25)-(3.28) 
and (3.30), summing over n = 0,1, 2,...,/ — 1 with 1 < I < N, and using the properties 
of the L^-projection and the regularity assumption (1.17), we obtain 


(3.31) 


1 

2 Ve 


-^+S^\T 


E 




+ 7® 
4 




n=l 


< -E 
- 2 


lie' 


0||2 


iL'^m 


+ -rE 


iive°-^111.(1,) 


+ —tE 
16 




C4 


+ C3(l +(52 +(5^) + ^ Tr 


C .2 


+ Co(1 + (52+(5^)T/i 2”“(’--*-1) + h|.r/j2min(r+l,s) 


^ + (5M E 


^ E liy Ili^CD) 

_ n=0 
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when ~ ^ ^.nd 


(3.32) 




E 


lie 


l \\2 

L 2 (D) 




■Eii^^' 


•n ||2 

LHV) 


n=l 


< 

“ 2 


lirili^d.) 


-rE 


2 I , ^^2 + C 4 


+ ( Cs (l + (^ +(^) + 


Tt 


16 


-rE 


liv^^lli.(^) 


+ (7o(l + )Th 


1 + - + E 


2 _|_ min (r,s —1) _|_ ^0_lL.^2.7^/}2 (r+l,s) 


1-1 

■Eiif 

n=0 


n II2 

LHV) 


Th^ 


when have not explicitly tracked the constant ||-^||ci(:d) 

or and some other constants in (3.31) and (|3.32[ ). 

Finally, estimates ( |3.12 ) and (3.14) follow from ( 3.3l{ -(3.32), the discrete Gron- 
wall’s inequality, the L^-projection properties, the fact = 0 and the triangle in¬ 
equality. □ 


Remark 6. Error estimates in Theorem 3.1 remain unehanged if we eonsider the 


modified seheme (3.10). In faet, we only need to eheek the fifth term on the right-hand 
side of d^rfl).- 


^2 r rt 

(3.33) E[T 5 ] = -yE 


- —E 
2 


^ ((divB - b) ■ {Vu{s) - Vu{tn+i)),C^Y ds 

((divB - b) ■ {Wri^+^ + e""^^) ds 


< Iixiiy^^yar^ + (1 + ||X||y(^))( 5 "E I Iie"+'lli.(i,) 


llx|iy(g)E I iivr+ii.(i,) 


"+16® 


livr+ii2(p) 


r. 


Henee the third term ^tE 


live^l 


2 

L 2 (D) 


on the right-hand side of (3.31) or (3.32) is 


replaeed by ^rE 


live^lp 


L 2 (I>) 


Remark 7. Error estimates (3.12) and (3.14) are optimal with respeet to r, and 
optimal in the H ^-no rm but suboptimal in the -norm with respeet to h. Erom the 
proof of Theorem 3A, the suboptimal estimate in the L‘^-norm is eaused by gradient- 


type noises, i.e., the existenee ofT^, T 5 and Tq on the right-hand side of (3.17). The 
proof in Theorem 3.1 relies on the p-th moment estimate driTi ) for the strong solu¬ 
tion. Convergenee rates are expeeted to be lower under weaker regularity assumptions. 
We also note that the error bounds depend exponentially on 1/e, whieh seems to be 
pessimistie. However, this is the ease even in the deterministie ease (i.e., S = Oj un¬ 
less the standard error estimate teehnique is replaeed by a mueh involved nonstandard 
teehnique as used in JEj/. We intend to address this issue in a future work. 

4. Numerical experiments. In this section we present a couple two-dimensional 
numerical experiments to gauge the performance of the proposed fully discrete finite 
element methods with r = 1, i.e., Vh is the linear finite element space. We also nu¬ 
merically examine the influence of noises on the dynamics of the numerical interfaces. 
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We consider the SPDE (1.10) on the square domain V = [—0.5, 0.5]^ with two 
different initial conditions, and in both tests we take X{x) = (p{x)[xi +X 2 , xi — ^ 2 ]^, 
where 


^{x) := 


_ 0.001 

0.092 _ |cc |2 


if \x\ < 0.3, 
if |x| > 0.3. 


For both tests, we take the Brownian motion step to be 1 x 10“^ and compute M = 500 
Monte Carlo realizations. In this section, we use a new notion to denote 

the numerical solution, where (5 > 0 controls the noise intensity, e > 0 is the diffuse 
interface width, h and r are the mesh sizes, and u; is a discrete sample. 


Next, we describe the algorithm which we use to solve the discrete problem (3.1). 
Let Nh = diml4 and be the nodal basis of I 4 . Denote by the coefficient 

vector corresponding to the discrete solution time = 


(n + l)r. Suppose ~ 

equivalent to 


(the other case is similar), (|3.1|) is then 


(4.1) 


M + r A 




,n+l 


Jn("' 


n+lN 


= 1 


(1 +J)Mu' 


2 


(5AWn+lC2U", 


where M,A are the standard mass and stiffness matrices, respectively. Ax is the 
weighted stiffness matrix whose (i,j) component is • X), N(u”+^) is 

the nonlinear term, (Ci)ij = ((divB — b) • V'0j,'0i), (C 2 )ij = (V'0j,'0i) and W is the 
discrete Brownian motion with increments AWn+i = Wn+i — Wn- Since we want 
to generate as many samples as possible in order to recover quantities of statistical 
interests, it will be expensive to use Newton’s method to solve (4.1). Instead, we 


employ a cheaper fixed point iteration, although it may converge slower than Newton’s 
method for a particular realization. However, if it is utilized well, more efficient overall 
solution algorithm can be designed when combined with the direct linear solver. 


Notice that the first coefficient matrix on the left-hand side of (4.1) is independent 
of n and cj, compute its Cholesky factorization 

M + t(^A+ = LL'. 

L will be stored and re-use for every Monte Carlo realization and at every time 
step. Only backward and forward substitutions are needed to execute the fixed point 
iteration, hence, resulting in considerable amount of computational saving. From this 
point of view, the algorithm based on the fixed point iteration outperforms that based 
on Newton’s method, especially for large sample size M. Indeed, we observed in our 
experiment that the fixed point iteration is much more efficient. 


Test 1. The initial condition uq is chosen as 


, . fd{x)\ 

uo [X) = tanh( 1, 

^ V 26 ^ 


where tanh(x) = (e^ —e ^)/(e^-he ^) and d{x) represents the signed distance between 
the point x and the ellipse 


0.04 0.01 
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In Table 


4.1 


we plot the expected values of the I/^([0,T], || • 11^,2 (d))- norm of the 
errors and rates of convergence of the time discretization for varying r with the fixed 
parameters (5 = 1 and e = 0.1. The numerical results confirm the theoretical result of 
Theorem 13.11 



Expected values of error 

Order of convergence 

r=0.008 

0.09895 


t=0.004 

0.06557 

0.5937 

r=0.002 

0.04472 

0.5521 

T=0.001 

0.03136 

0.5120 


Table 4.1 

Computed time diseretization errors and eonvergenee rates 


In Figure [TTf|4.3[ we display a few snapshots of the zero-level set of the expected 
value of the numerical solution 




M 

\ ^ .,6,e,h,r 


M ^ 

i=l 




at several time points with e = 0.01, and three different noise intensity parameters 
6 = 0.1,1,10. We observe that the shape of the zero-level set of the expected value 
of the numerical solution undergoes more changes as S increases. 




0.4 - 

0.3 . 

0.2 - 

0.1 - 

0.4 

0.3 

0.2 

0.1 


o 

0 

o 

0.1 - 

0.2 . 

0.3 - 

0.4 - 

0!? _1_ 

-0.1 

-0.2 

-0.3 

-0.4 

-n ^ 



Fig. 4.1. Snapshots of the zero-level set of at time t = 0, 0.020, 0.040, 0.043 with 6 = 0.1 

and e = 0.01. 

Next, we study the effect of e on the evolution of the numerical interfaces. To 
the end, we fix (5 = 0.1. In Figure |T^ we depict four snapshots at four fixed time 
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Fig. 4.2. Snapshots of the zero-level set of at time t = 0.005,0.020,0.040,0.043 with 

5 = 1 and e = 0.01. 

points of the zero-level set of the expected value of the numerical solution at 

three different e = 0.01,0.011,0.02. The numerical interface clearly converges to the 
stochastic mean curvature flow as e ^ 0, and it evolves faster in time for larger e. 
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Fig. 4.3. Snapshots of the zero-level set of at time t 

0.0025, 0.0050, 0.0100, 0.0200, 0.0250, 0.0280 with S = 10 and e = 0.01. 


Test 2. First, we define 


tanh(^) := 

i’2{x2) : = 

i’4{x2) : = 


py — p y 


ey + e-y 

1 - V1.92X2 +0.2304 
2 ’ 
1-V-1-92X2 +0.2304 


i^l{x2) 

i^3{x2) 

'05(^2) 


— 1 + \/0i8x2^^^i04 

2 ’ 

— 1 + V—0.8x2 + 0.04 


1 - 0.2451x^ 
0.0049 


and consider the initial condition uo{xi^X 2 ) = iii(3xi, 3x2), where ixi(xi, X 2 ) is defined 
by (cf. Test 2 in [9]): 
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Fig. 4.4. Snapshots of the zero-level set of at time t = 0.0010,0.0050,0.0125,0.0175 

with S = 0.1 and e = 0.01, 0.011, 0.02. 


Ui{xi,X2) 


tanh(^(-y(xi - 0.14)2 - 0.15)2)), 

tanh(^(-^(xi - 0.14)2 + 0.15)2)), 

tanh(^(-A/(xi + 0.3)2 - 0.15)2)), 

tanh(^(-^(xi + 0.3)2 O.IS)^)), 

tanh(^(y(a;i - 0.5)2 ^2 _ 0.39)), 

tanh(^(y^(xi + 0.5)2 _|_ ^2 _ o.25)), 

-0.15)), 


tanh(^(^(xi - 0.5)2 ^2 _ 0.39)), 

tanh(^(y^(xi - 0.5)^ + - 0.39)), 

tanh(^(^(xi + 0.5)2 ^2 _ o.25)), 

tanh(^(y(a;i + 0.5)2 ^ ^2 _ o.25)). 


if xi > 0.14,0 <X2< —^{xi — 0.5), 
if xi > 0.14, ^(xi — 0.5) < X2 < 0, 
if xi < —0.3, {) < X2 < \{xi -\- 0.5), 
if xi < —0.3, —|(xi + 0.5) < X2 < 0, 
if xi > 0.14,0:2 > —^{xi — 0.5) 
or X 2 < ^{xi — 0.5), 
if xi < —0.3, X2 > —|(xi + 0.5) 
or 0:2 < -\{xi + 0.5), 
if -0.3 < xi < 0.14, 

^1{X2) <Xi< I/^2(X2) 

and '03(^r2) ^ xi < '04(a:2), 
if — 0.3 < xi < 0.14, xi > '02(iC2) 
and xi > '05(x2), 

if — 0.3 <xi< 0.14, xi > '04(iC2) 
and xi > '05(^r2), 

if — 0.3 <xi< 0.14, xi < '01 (X2) 
and xi < 05(0:2), 

if — 0.3 < 0:1 < 0.14, 0:1 < 03(iC2) 
and 0:1 < 05(3:2). 
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We note that the initial condition is not smooth due to the dumbbell shape of 
the zero-level set. Nevertheless, we study the effects of S and e on the evolution of 
numerical interfaces. Figure [4^4 . 7| display a few snapshots of the zero-level set of the 
expected value of the numerical solution at several time points with e = 0.01 

and three different noise intensity 5 = 0.1,1,10. Similar to Test 1, the zero-level set 
evolves faster and the shape undergoes more changes for larger 5. Figure |4^ plots 
snapshots at four fixed time point of the zero-level set of with = 0.1 and 

e = 0.01,0.011,0.02. Again, it suggests the convergence of the zero-level set at each 
time point to the stochastic mean curvature flow as e ^ 0. 






Fig. 4.5. Snapshots of the zero-level set of at time t = 0, 0.040, 0.200, 0.456 with 6 = 0.1 

and e = 0.01. 
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Fig. 4.7. Snapshots of the zero-level set of at time t 

0.004, 0.040, 0.080, 0.140, 0.180, 0.216 with S = 10 and e = 0.01. 
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Fig. 4.8. Snapshots of the zero-level set of at time t = 0.008,0.040,0.080,0.120. 
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